Heterogeneity in intrahepatic macrophage populations and druggable target expression in patients with steatotic liver disease-related fibrosis.

Background & Aims: Clinical trials for reducing fibrosis in steatotic liver disease (SLD) have targeted macrophages with variable results. We evaluated intrahepatic macrophages in patients with SLD to determine if activity scores or fibrosis stages influenced phenotypes and expression of druggable targets, such as CCR2 and galectin-3. Methods: Liver biopsies from controls or patients with minimal or advanced fibrosis were subject to gene expression analysis using nCounter to determine differences in macrophage-related genes (n = 30). To investigate variability among individual patients, we compared additional biopsies by staining them with multiplex antibody panels (CD68/CD14/CD16/CD163/Mac387 or CD163/CCR2/galectin-3/Mac387) followed by spectral imaging and spatial analysis. Algorithms that utilize deep learning/artificial intelligence were applied to create cell cluster plots, phenotype profile maps, and to determine levels of protein expression (n = 34). Results: Several genes known to be pro-fibrotic (e.g. CD206, TREM2, CD163, and ARG1) showed either no significant differences or significantly decreased with advanced fibrosis. Although marked variability in gene expression was observed in individual patients with cirrhosis, several druggable targets and their ligands (e.g. CCR2, CCR5, CCL2, CCL5, and LGALS3) were significantly increased when compared to patients with minimal fibrosis. Antibody panels identified populations that were significantly increased (e.g. Mac387+), decreased (e.g. CD14+), or enriched (e.g. interactions of Mac387) in patients that had progression of disease or advanced fibrosis. Despite heterogeneity in patients with SLD, several macrophage phenotypes and druggable targets showed a positive correlation with increasing NAFLD activity scores and fibrosis stages. Conclusions: Patients with SLD have markedly varied macrophage- and druggable target-related gene and protein expression in their livers. Several patients had relatively high expression, while others were like controls. Overall, patients with more advanced disease had significantly higher expression of CCR2 and galectin-3 at both the gene and protein levels. Impact and implications: Appreciating individual differences within the hepatic microenvironment of patients with SLD may be paramount to developing effective treatments. These results may explain why such a small percentage of patients have responded to macrophage-targeting therapies and provide additional support for precision medicine-guided treatment of chronic liver diseases.

Patients with fibrosis from non-alcoholic steatohepatitis have heterogeneous intrahepatic macrophages and therapeutic targets.

Background and Aims: In clinical trials for reducing fibrosis in NASH patients, therapeutics that target macrophages have had variable results. We evaluated intrahepatic macrophages in patients with non-alcoholic steatohepatitis to determine if fibrosis influenced phenotypes and expression of CCR2 and Galectin-3. Approach & Results: We used nCounter to analyze liver biopsies from well-matched patients with minimal (n=12) or advanced (n=12) fibrosis to determine which macrophage-related genes would be significantly different. Known therapy targets (e.g., CCR2 and Galectin-3) were significantly increased in patients with cirrhosis.However, several genes (e.g., CD68, CD16, and CD14) did not show significant differences, and CD163, a marker of pro-fibrotic macrophages was significantly decreased with cirrhosis. Next, we analyzed patients with minimal (n=6) or advanced fibrosis (n=5) using approaches that preserved hepatic architecture by multiplex-staining with anti-CD68, Mac387, CD163, CD14, and CD16. Spectral data were analyzed using deep learning/artificial intelligence to determine percentages and spatial relationships. This approach showed patients with advanced fibrosis had increased CD68+, CD16+, Mac387+, CD163+, and CD16+CD163+ populations. Interaction of CD68+ and Mac387+ populations was significantly increased in patients with cirrhosis and enrichment of these same phenotypes in individuals with minimal fibrosis correlated with poor outcomes. Evaluation of a final set of patients (n=4) also showed heterogenous expression of CD163, CCR2, Galectin-3, and Mac387, and significant differences were not dependent on fibrosis stage or NAFLD activity. Conclusions: Approaches that leave hepatic architecture intact, like multispectral imaging, may be paramount to developing effective treatments for NASH. In addition, understanding individual differences in patients may be required for optimal responses to macrophage-targeting therapies.

Subclonal evolution and expansion of spatially distinct THY1-positive cells is associated with recurrence in glioblastoma.

PURPOSE: Glioblastoma(GBM) is a lethal disease characterized by inevitable recurrence. Here we investigate the molecular pathways mediating resistance, with the goal of identifying novel therapeutic opportunities. EXPERIMENTAL DESIGN: We developed a longitudinal in vivo recurrence model utilizing patient-derived explants to produce paired specimens(pre- and post-recurrence) following temozolomide(TMZ) and radiation(IR). These specimens were evaluated for treatment response and to identify gene expression pathways driving treatment resistance. Findings were clinically validated using spatial transcriptomics of human GBMs. RESULTS: These studies reveal in replicate cohorts, a gene expression profile characterized by upregulation of mesenchymal and stem-like genes at recurrence. Analyses of clinical databases revealed significant association of this transcriptional profile with worse overall survival and upregulation at recurrence. Notably, gene expression analyses identified upregulation of TGFß signaling, and more than one-hundred-fold increase in THY1 levels at recurrence. Furthermore, THY1-positive cells represented <10% of cells in treatment-naïve tumors, compared to 75-96% in recurrent tumors. We then isolated THY1-positive cells from treatment-naïve patient samples and determined that they were inherently resistant to chemoradiation in orthotopic models. Additionally, using image-guided biopsies from treatment-naïve human GBM, we conducted spatial transcriptomic analyses. This revealed rare THY1+ regions characterized by mesenchymal/stem-like gene expression, analogous to our recurrent mouse model, which co-localized with macrophages within the perivascular niche. We then inhibited TGFBRI activity in vivo which decreased mesenchymal/stem-like protein levels, including THY1, and restored sensitivity to TMZ/IR in recurrent tumors. CONCLUSIONS: These findings reveal that GBM recurrence may result from tumor repopulation by pre-existing, therapy-resistant, THY1-positive, mesenchymal cells within the perivascular niche.

MHC Class II Expression Influences the Composition and Distribution of Immune Cells in the Metastatic Colorectal Cancer Microenvironment.

Despite advances in therapy over the past decades, metastatic colorectal cancer (mCRC) remains a highly morbid disease. While the impact of MHC-I on immune infiltration in mCRC has been well studied, data on the consequences of MHC-II loss are lacking. Multiplex fluorescent immunohistochemistry (mfIHC) was performed on 149 patients undergoing curative intent resection for mCRC and stratified into high and low human leukocyte antigen isotype DR (HLA-DR) expressing tumors. Intratumoral HLA-DR expression was found in stromal bands, and its expression level was associated with different infiltrating immune cell makeup and distribution. Low HLA-DR expression was associated with increased intercellular distances and decreased population mixing of T helper cells and antigen-presenting cells (APC), suggestive of decreased interactions. This was associated with less co-localization of tumor cells and cytotoxic T lymphocytes (CTLs), which tended to be in a less activated state as determined by Ki67 and granzyme B expression. These findings suggest that low HLA-DR in the tumor microenvironment of mCRC may reflect a state of poor helper T-cell interactions with APCs and CTL-mediated anti-tumor activity. Efforts to restore/enhance MHC-II presentation may be a useful strategy to enhance checkpoint inhibition therapy in the future.

Stress response gene family expansions correlate with invasive potential in teleost fish.

The bluegill sunfish Lepomis macrochirus and the closely related redear sunfish Lepomis microlophus have important ecological and recreational value and are widely used for research and aquaculture. While both species have been introduced outside of their native ranges, only the bluegill is considered invasive. Here, we report de novo transcriptome assemblies for these fish as a resource for sunfish biology. Comparative analyses of the transcriptomes revealed an unexpected, bluegill-specific expansion in the HSP70 and HSP90 molecular chaperone gene families. These expansions were not unique to the bluegill as expansions in HSP70s and HSP90s were identified in the genomes of other teleost fish using the NCBI RefSeq database. To determine whether gene family expansions are specific for thermal stress responses, GST and SOD gene families that are associated with oxidative stress responses were also analyzed. Species-specific expansions were also observed for these gene families in distinct fish species. Validating our approach, previously described expansions in the MHC gene family were also identified. Intriguingly, the number of HSP70 paralogs was positively correlated with thermotolerance range for each species, suggesting that these expansions can impact organismal physiology. Furthermore, fish that are considered invasive contained a higher average number of HSP70 paralogs than non-invasive fish. Invasive fish also had higher average numbers of HSP90, MHC and GST paralogs, but not SOD paralogs. Taken together, we propose that expansions in key cellular stress response gene families represent novel genetic signatures that correlate with invasive potential.

Developing and deploying an integrated workshop curriculum teaching computational skills for reproducible research.

Inspired by well-established material and pedagogy provided by The Carpentries (Wilson, 2016), we developed a two-day workshop curriculum that teaches introductory R programming for managing, analyzing, plotting and reporting data using packages from the tidyverse (Wickham et al., 2019), the Unix shell, version control with git, and GitHub. While the official Software Carpentry curriculum is comprehensive, we found that it contains too much content for a two-day workshop. We also felt that the independent nature of the lessons left learners confused about how to integrate the newly acquired programming skills in their own work. Thus, we developed a new curriculum that aims to teach novices how to implement reproducible research principles in their own data analysis. The curriculum integrates live coding lessons with individual-level and group-based practice exercises, and also serves as a succinct resource that learners can reference both during and after the workshop. Moreover, it lowers the entry barrier for new instructors as they do not have to develop their own teaching materials or sift through extensive content. We developed this curriculum during a two-day sprint, successfully used it to host a two-day virtual workshop with almost 40 participants, and updated the material based on instructor and learner feedback. We hope that our new curriculum will prove useful to future instructors interested in teaching workshops with similar learning objectives.

CGAT: Cell Graph ATtention Network for Grading of Pancreatic Disease Histology Images.

Early detection of Pancreatic Ductal Adenocarcinoma (PDAC), one of the most aggressive malignancies of the pancreas, is crucial to avoid metastatic spread to other body regions. Detection of pancreatic cancer is typically carried out by assessing the distribution and arrangement of tumor and immune cells in histology images. This is further complicated due to morphological similarities with chronic pancreatitis (CP), and the co-occurrence of precursor lesions in the same tissue. Most of the current automated methods for grading pancreatic cancers rely on extensive feature engineering involving accurate identification of cell features or utilising single number spatially informed indices for grading purposes. Moreover, sophisticated methods involving black-box approaches, such as neural networks, do not offer insights into the model's ability to accurately identify the correct disease grade. In this paper, we develop a novel cell-graph based Cell-Graph Attention (CGAT) network for the precise classification of pancreatic cancer and its precursors from multiplexed immunofluorescence histology images into the six different types of pancreatic diseases. The issue of class imbalance is addressed through bootstrapping multiple CGAT-nets, while the self-attention mechanism facilitates visualization of cell-cell features that are likely responsible for the predictive capabilities of the model. It is also shown that the model significantly outperforms the decision tree classifiers built using spatially informed metric, such as the Morisita-Horn (MH) indices.

Teaching Python for Data Science: Collaborative development of a modular & interactive curriculum.

We are bioinformatics trainees at the University of Michigan who started a local chapter of Girls Who Code to provide a fun and supportive environment for high school women to learn the power of coding. Our goal was to cover basic coding topics and data science concepts through live coding and hands-on practice. However, we could not find a resource that exactly met our needs. Therefore, over the past three years, we have developed a curriculum and instructional format using Jupyter notebooks to effectively teach introductory Python for data science. This method, inspired by The Carpentries organization, uses bite-sized lessons followed by independent practice time to reinforce coding concepts, and culminates in a data science capstone project using real-world data. We believe our open curriculum is a valuable resource to the wider education community and hope that educators will use and improve our lessons, practice problems, and teaching best practices. Anyone can contribute to our Open Educational Resources on GitHub.

GWAS of thyroid stimulating hormone highlights pleiotropic effects and inverse association with thyroid cancer.

Thyroid stimulating hormone (TSH) is critical for normal development and metabolism. To better understand the genetic contribution to TSH levels, we conduct a GWAS meta-analysis at 22.4 million genetic markers in up to 119,715 individuals and identify 74 genome-wide significant loci for TSH, of which 28 are previously unreported. Functional experiments show that the thyroglobulin protein-altering variants P118L and G67S impact thyroglobulin secretion. Phenome-wide association analysis in the UK Biobank demonstrates the pleiotropic effects of TSH-associated variants and a polygenic score for higher TSH levels is associated with a reduced risk of thyroid cancer in the UK Biobank and three other independent studies. Two-sample Mendelian randomization using TSH index variants as instrumental variables suggests a protective effect of higher TSH levels (indicating lower thyroid function) on risk of thyroid cancer and goiter. Our findings highlight the pleiotropic effects of TSH-associated variants on thyroid function and growth of malignant and benign thyroid tumors.

Mathematical Modeling of the Metastatic Colorectal Cancer Microenvironment Defines the Importance of Cytotoxic Lymphocyte Infiltration and Presence of PD-L1 on Antigen Presenting Cells.

BACKGROUND: Although immune-based therapy has proven efficacious for some patients with microsatellite instability (MSI) colon cancers, a majority of patients receive limited benefit. Conversely, select patients with microsatellite stable (MSS) tumors respond to checkpoint blockade, necessitating novel ways to study the immune tumor microenvironment (TME). We used phenotypic and spatial data from infiltrating immune and tumor cells to model cellular mixing to predict disease specific outcomes in patients with colorectal liver metastases. METHODS: Formalin fixed paraffin embedded metastatic colon cancer tissue from 195 patients were subjected to multiplex immunohistochemistry (mfIHC). After phenotyping, the G-function was calculated for each patient and cell type. Data was correlated with clinical outcomes and survival. RESULTS: High tumor cell to cytotoxic T lymphocyte (TC-CTL) mixing was associated with both a pro-inflammatory and immunosuppressive TME characterized by increased CTL infiltration and PD-L1+ expression, respectively. Presence and engagement of antigen presenting cells (APC) and helper T cells (Th) were associated with greater TC-CTL mixing and improved 5-year disease specific survival compared to patients with a low degree of mixing (42% vs. 16%, p = 0.0275). Comparison of measured mixing to a calculated theoretical random mixing revealed that PD-L1 expression on APCs resulted in an environment where CTLs were non-randomly less associated with TCs, highlighting their biologic significance. CONCLUSION: Evaluation of immune interactions within the TME of metastatic colon cancer using mfIHC in combination with mathematical modeling characterized cellular mixing of TCs and CTLs, providing a novel strategy to better predict clinical outcomes while identifying potential candidates for immune based therapies.